Cyclic bases 
of zero-curvature representations: 
five illustrations to one concept 

S. Yu. Sakovich 

Institute of Physics, National Academy of Sciences, 220072 Minsk, Belarus 

e-mail: saks@pisem.net 

Abstract 

The paper contains five examples of using cyclic bases of zero- 
curvature representations in studies of weak and strong Lax pairs, 
hierarchies of evolution systems, and recursion operators. 



1 Introduction 

Many remarkable equations of nonlinear mathematical physics can be repre- 
sented as the compatibility condition 

Z = D t X - D X T + [X,T] = (1) 

of the overdetermined linear system 

* a = XV, % = Ttf , (2) 

where D x and D t stand for the total derivatives; X and T are n x n matrix 
functions of the independent variables x and t, dependent variables u % (x, t) 
and derivatives of u 1 ; ^ (x, t) is an n-component column; the square brackets 
denote the matrix commutator. The condition (|I|), usually referred to as the 
zero-curvature representation (ZCR), is said to represent a given equation if 
any solution u % of the equation satisfies ([!]). The linear transformation 

V(x,t) =G*(x,t) (3) 
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of the auxiliary vector function \l/ with any nondegenerate n x n matrix 
function G (x, t, u\ u l x , . . .), det G ^ 0, generates the gauge transformation of 
the matrices X and T: 

X' = GXG- 1 + (D X G) G-\ T' = GTG- 1 + (D t G) G~ l . (4) 

Since the matrix Z in (H) is transformed under as 

Z' = GZG~\ (5) 

any two ZCRs related by a gauge transformation (^j) are considered as equiv- 
alent (see, e.g., JTJ]). The existence of such an equivalence between ZCRs 
suggests that ZCRs should be studied by gauge-invariant methods. 

In the present paper, we give five examples illustrating the use of the 
cyclic bases of ZCRs — one of the central concepts of our work 0, where a 
gauge-invariant description of ZCRs of evolution equations was proposed^. 
The paper is organized as follows. In Section |2|, we study a continual class 
of second-order evolution equations admitting weak Lax pairs, and construct 
an extension of this class to any order. In Section ^, we consider removability 
of a parameter from a weak Lax pair of the Burgers equation, and then, in 
Section [|, we derive the Burgers equation's recursion operator from the cyclic 
basis of a strong Lax pair. In Section |j, we apply this method of deriving 
recursion operators from ZCRs to a system of coupled KdV equations, in 
order to show how the method works with systems, and then, in Section |6|, 
we use the same method again, in order to explain a strange structure of a 
recursion operator of one KdV-mKdV system. Section [F] contains concluding 
remarks. 



2 Second-order evolution equations 

Recently, Marvan |j| proved that every second-order scalar evolution equa- 
tion possessing an irreducible sl 2 -valued ZCR can be brought by a contact 
transformation into the form 

u t = p x u 2 u xx + 2f3 xx u 2 u x + Af3u x + {f3 xxx - Af3 x ) u 3 - Af3 x u, (6) 

1 A much more general (but very abstract, indeed) gauge-invariant description, appli- 
cable to ZCRs of any non-overdetermined systems of PDEs, was introduced by Marvan in 
later, in Marvan gave explicit expressions for the objects involved. 
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where j3 (x, t) is any function with (3 X ^ 0, and then the ZCR of (§) is given 
by (P with 



X= u (7) 
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One may wonder why the class of evolution equations (§) contains the arbi- 
trary function (3 (x, t). The origin of j3 (x, t) can be revealed by the technique 
of cyclic bases of ZCRs @. 

Let us solve the problem of finding all the (1 + l)-dimensional scalar 
evolution equations 



l^t / tj ^xi ■ ■ ■ i 'U'x...x) 



(9) 



which admit ZCRs (|T]) with the predetermined matrix X ([?D and any 2x2 
matrices T (x, t, u, u x , . . . , u x _ x ) (traceless, without loss of generality), the 
highest orders of derivatives u x ,„ x m / and T being not fixed. Since X 
(0) contains no parameter, this problem is equivalent to finding the complete 
class of evolution equations possessing weak Lax pairs with the predetermined 
spatial part. 

First, we rewrite the ZCR of (|D|) in its equivalent (characteristic) form 



fC = V X T, 



(10) 



where C is the characteristic matrix of the ZCR (in the present case, C is 
simply C = dX/du, because X (^) contains no derivatives of u @, ||), and 
the operator V x is defined by V X Y = D X Y — [X, Y] for any (here, 2x2) 
matrix Y. 

Second, we compute C, V X C, V X C and V X C for the matrix X (^), 
and find that the cyclic basis (i.e. the maximal sequence of linearly inde- 
pendent matrices V X C, i = 0,1,2,...) is three-dimensional in this case, 
{C, V X C, V X C}, and that the closure equation for the cyclic basis is 



VlC = a x C + a 2 V x C + a 3 V 2 x C 



'ir, 
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with 

8u x — 2v? x 8u x u xx V2u x — 2u xxx 
Cb\ = ; 5 1 , 

4 — lCh4 6u xx 7u x . . 

a 2 = ^ - + 4, a 3 = x -. 12 

xir u u 

Third, we decompose the unknown matrix T over the cyclic basis as 

T = b 1 C + b 2 V x C + b 3 V 2 x C, (13) 

where b\, b 2 and b 3 are unknown scalar functions of x, t, u, u x , . . . , u x ___ x . Sub- 
stituting (|IB|) into ( PH| ) and using the closure equation ([□]) and the fact of 
linear independence of C, V X C and V^C, we obtain the following relations: 

b 2 = -D x b 3 - a 3 b 3 , b x = -D x b 2 - a 2 b 3 , f = D x bi + a x b 3 , (14) 

where the function 63 remains arbitrary. 

Fourth, we combine (0), (|13|) and (|I4D, use the explicit form of C, V ' X C 
and V^C, and thus obtain the following expressions for T and / in terms of 
one arbitrary function b 3 = p (x, t, u, u x , . . . , u x ,„ x ), where the order of u x ,„ x 
is also arbitrary: 



with 



, 1 ^9 5u x „ 4 — 9ul + 3uu Xx , 
T u = --Dl + -^D x + x — x - p, 



-T 12 = ^rAr h s — p, 



9 * 1 ^ 



= \-1d. + 1±P±) I >. (16) 

IT It 



and 



f - / n3 _ J^E n2 _ 4 + 4w 2 - 24m 2 + 8uu xx 

It IT 



lQu x + \2u z u x — 36u^ + 27uu x u xx — 3u 2 u xxx 
w 



^.^^ x 1 — ^ v^j; ^wv^„ I HI UjU> X U> XX ^ XXX \ 
+ -3 IP' ( 17 ) 



The problem has been solved. Every evolution equation (0), which be- 
longs to the continual class determined by ( |TT| ) with any function p(x, t, u 



u x , . . . , u x ,,, x ) of any order in u x ,„ x , admits the ZCR ([!]) with X given by (|7|) 
and T determined by (|15|)-(|16|). As it was shown in ||, such a continual 
class of evolution equations, which all possess ZCRs with a predetermined 
matrix X, exists for every matrix X containing no parameter. Of course, the 
studied case of X (0) is not an exception. 

Now we can see that the origin of the arbitrary function f3 (x, t) in (|6|) is 
the arbitrary function p (x,t,u,u x , . . . ,u x ,„ x ) in (|17|) . Indeed, if we require 
that / ([T7|) is a function of x, t, u, u x and u xx only, i.e. that the order of 
the represented evolution equations is two, then p = p(x, t, u) (otherwise / 
contains u x ___ x of order four or higher) and dp/du = 3p/u (this follows from 
df /du xxx = 0), that is 

p = P(x,t)u 3 (18) 

with any function j3 {(3 X ^ if we need df /du xx ^ 0). Finally, the equation 
(HD relates (0) with (|), as well as ©-(H) with (|). 



3 Burgers equation's weak Lax pair 

Our next example is the ZCR ([!]) with 



, \r\u \u x + \u 2 + \r\u i 

T = L \ ' x 4 \ 4 , (20) 

where rj is a parameter. This ZCR of the Burgers equation u t = u xx + uu x 
was introduced by Cavalcante and Tenenblat 0. The parameter rj, however, 
is not an essential ('spectral') parameter: it can be removed ('gauged out') 
from (plj) and (|20|) by a gauge transformation (|4]) (therefore this ZCR is 



a weak Lax pair of the Burgers equation). More precisely, it was recently 
shown by Marvan [7] that any nonzero value of 77 in these X and T can be 
changed into any other nonzero value (e.g., into rj = 1) by an appropriate 
gauge transformation.^ 

Let us see how the technique of cyclic bases allows to discover, in a short 
computational way, that the case of 77 = is essentially different from any 

2 The fact of removability of 77 from X (|l9|) was also noted in j2|, in slightly different 
notations, in relation with ZCRs of third-order evolution equations B. 
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other case of nonzero 77, and that the parameter 77 should be expected to 
be removable from this ZCR. The necessary information can be obtained 
from the dimension of the cyclic basis and from the coefficients of the closure 
equation. 

Following the scheme of Section we find the characteristic matrix and 
its lower-order covariant derivatives to be 

V^=2(-' -/), V^f(-/ -/). (2!) 

We see from ( PT| ) that, when rj ^ 0, the cyclic basis is two-dimensional, 
{C, V X C}, with the closure equation 

11 

V 2 X C = -V.C, (22) 

whereas in the case of rj = the cyclic basis is one-dimensional, consisting of 
the characteristic matrix C itself, with the closure equation 

V X C = 0. (23) 

Since the matrices V^C (i = 0, 1, 2, . . .) are transformed under (|j) as 
V^C" = G{V x C)G~ l |@, U, the dimension of the cyclic basis and the co- 
efficients of the closure equation are invariants with respect to gauge trans- 
formations. Therefore no gauge transformation can change the dimension 
of the cyclic basis and relate the case of rj ^ with the case of 77 = 0. 
Moreover, the coefficients of the closure equation (|22|) do not contain 77, and 
this is highly indicative of the parameter removability^. On the contrary, if 
there is a parameter in a coefficient of a closure equation of a ZCR, then the 
parameter is essential, not removable from the ZCR, for the reason of the 
gauge invariance of coefficients of closure equations, and Sections £|-§] provide 
examples illustrating this phenomenon. 

Using the method described in Section Q, we can construct continual 
classes of evolution equations admitting ZCRs with X given by (|19|) , sepa- 
rately for the two cases of the cyclic basis structure: 

u t = D xP for 77 = 0, (24) 



3 An easy way to prove this removability of 77 is to solve (Q) with respect to G directly, 
taking X' = X\ n= i and T' = T| ??= i and starting from the evident necessary conditions 
G = G (x, t) and [C, G] = 0. A different approach was used in 0. 
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Ut = D x (D x + ^p for 77 # 0, (25) 

where p is any function of x,t,u,u x , . . . ,u x ___ x , of any order in u x _ x . Note 
that the orders of the linear differential operators in the right-hand sides of 
(|24| ) and fl2"5]) are equal to the dimensions of the corresponding cyclic bases, 
in accordance with the general theorem on continual classes ||. 

The class ( ^ ) consists of evolution equations admitting the conserved 
density u. This is a natural consequence of the cyclic basis structure in 
the case of rj = 0. Indeed, if we consider a (1 x l)-dimensional ZCR with 
X = u and any function T, then this ZCR is simply a conservation law, 
D t X — D X T = 0, due to [X, T] = 0, and the closure equation for X = u is 
exactly (p3|). 

The class ( P5|) is more narrow than (^4j) but much wider than habitual 
discrete hierarchies of integrable equations. See pi for the results of the 
singularity analysis which suggest that most of equations in such continual 
classes are non-integrable. 

The fact that the Burgers equation u t = u xx + uu x admits a ZCR with 
X QT9D is equivalent to the fact that the equation belongs to the classes 
and (p5|). It tells nothing about the Burgers equation's integrability. 



4 Burgers equation's recursion operator 

Now, let us study a strong Lax pair of the Burgers equation 

u t = u xx + 2uu x , (26) 
namely, the ZCR (|I|) with 

x -{ i -Uu + \)h (27) 



T=(^ + U \-^ ' \ A, (28) 
V u-X [u x + u — X ) J v ' 

where A is any constant. ^ Note that these X and T are lower triangular ma- 
trices. However, contrary to the approach adopted in 0, we do not exclude 
such lower triangular (or reducible) ZCRs from our consideration. 

4 We are unable to indicate the first appearance of this ZCR in the literature. In implicit 
form, this ZCR can be found in ||, among non-Abelian pseudopotentials of the Burgers 
equation. 
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Computing C, V ' X C and V 2 X C for the matrix X (p7|), we find that the 
cyclic basis is two-dimensional, {C,V X C}, with the closure equation 

V 2 C = (u + A) V^C. (29) 

The presence of A in the gauge-invariant coefficient u + A of (|29|) proves that 
the parameter A is essential, not removable by gauge transformations. The 
matrix T (ESI) can be decomposed over the cyclic basis as 



T = (u x + u 2 - A 2 ) C + (A - u) V X .C. (30) 

Let us find all the evolution equations (|D|) which admit ZCRs ([!]) with X 
determined by ( p7|) and T decomposable over the cyclic basis as 



T = b 1 C + b 2 V x C, (31) 

where hi and 6 2 are any scalar functions of x, t, u, u x , . . . , and A. Using 
fllOl) , (|3~ll), (|29| ) and the fact of linear independence of C and V X C, we find 
that 

b 1 = (D x + u + X)p, (32) 



/ = (D 2 X + (u + A) A, + p, (33) 

where p (A, x, t, u, u x , . . . , u x „, x ) = —b<i is any function, of any order in u x „, x . 
Note that the linear differential operator in the right-hand side of (^) con- 
tains the parameter A. If A is any fixed constant, then our problem is com- 
pletely solved by the expression ( |3"B| ) for the right-hand side of every repre- 
sented equation @, and we obtain a typical continual class. However, A is 
a free parameter, and we have to determine the set of functions p satisfying 
the condition df /d\ = imposed on (^). Here we give a detailed derivation 
of those admissible p, in order to omit details of the same way of reasoning 
in Sections |5] and [| 

First, we introduce the operators 

M = D 2 X + uD x + u X} N = D x , (34) 

expand the function p as 

p = po + Api + X 2 p 2 + X 3 p 3 + • • • , (35) 
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where pi (i = 0, 1, 2, . . .) are unknown functions of x, t, u, u x , . . . , u x ,,, x , and 
thus rewrite ([33]) in the formf] 

/ = (M + XN) (po + Xpi + X 2 P 2 + •••)• ( 36 ) 



Second, we obtain from ( |5B| ) the expression for / satisfying the condition 
df/d\ = 0, 

f = Mp , (37) 
as well as the recursion relation determining the functions pi, 

Mp i+1 + Npi = 0, i = 0,1,2, .... (38) 



Third, taking into account that M and N ( p4|) are differential operators 
of orders two and one, respectively, we conclude from (|3~8"D that there exists a 
sufficiently large number m, such that p m does not contain « and derivatives 
of u. Then the recursion relation (j38"|) for i > m and the explicit expressions 
for M and iV fl3"I|) lead to p m = <fi (t) , with any 0, and = for i > m. 

The problem has been solved: the right-hand side of every evolution 
equation sought is determined by fl3"T|), where po is given by the recursion 
( |38| ) starting from p m = (t), with any function <fi (t) and any integer m > 0, 
the operators M and N being defined by fl34j) . 

However, fourth, the result can be represented in a more compact form. 
Using the operator iV _1 , formally inverse^ to N, we express po from (^) as 
Po = (-N- l M) m p m , rewrite (0) as / = M (— A^ _1 M) m N~ l Np m , take into 
account that iVp m = 0, and thus obtain 

/ = .R^+iO, (39) 

where 

R = MN- 1 =D x + u + u x D~ l . (40) 



Consequently, the matrix X (|27|) is related not to a continual class of 
evolution equations, as it would be if there was no essential parameter in 
X, but to the discrete hierarchy determined by (0), which consists of the 
evolution equations 

Ut = 01 (t) U x + 02 (t) (U xx + 2UU X ) 

+ 03 (t) (u xxx + 3uu xx + 3u 2 x + 3u 2 u x ) H (41) 



5 If a singularity is suspected in p at A = 0, one may use ( p6[ ) after an infinitesimal shift 
of A, A i-> A + e, M M - eN, N i-> AT, e 0, with no effect on the final result (|o|). 

6 In the sense that N~ 1 a denotes any function b (x, t, u, u x , . . . , u x ... x ) such that a = Nb 
if b exists. 
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with any functions 0i, 02, etc. All the equations (^TJ) can be linearized by the 
Cole-Hopf transformation. The obtained recursion operator R (j^), which 
generates the hierarchy of equations fl4T|), is well known as a recursion opera- 
tor of the Burgers equation in the sense of higher symmetries (or generalized 
symmetries, or Lie-Backlund algebras) (see, e.g., flOU). 



5 Coupled KdV equations 

The method, used in Section |] to derive the Burgers equation's hierarchy 
and recursion operator, can be applied to systems of evolution equations as 
well. Therefore our next example is 

Ut = u xxx — 12uu x , v t = 4:V XXX — 6vu x — 12uv x . (42) 



as a new 



These coupled KdV equations were found by Giirses and Karasu |TTJ 
system possessing a second-order recursion operator (in the sense of higher 
symmetries). Later, the system (|42|) appeared in |f[^| , in a list of coupled 
KdV equations that passed the Painleve test for integrability well, and its 
ZCR (HD with 3x3 matrices X and T containing an essential parameter was 
also found there. For what follows, we need to know the explicit form of the 



matrix X from JT^j . 




X = 2 , (43) 



where a is a parameter. No explicit expression for the matrix T is necessary: 
it is sufficient to know that T can be decomposed over the cyclic basis. 

In this case of a system of two equations, we have two characteristic ma- 
trices, C u and C v , which are simply C u = dX/du and C v = dX/dv because 
X ( |4~3|) contains no derivatives of u and v p], |J. Computing V X C U , V X C V , 
V 2 X C U , V 2 X C V and V X C U , we find that the cyclic basis is five-dimensional, 
{C u , C v , V X C U , V X C V , V^C U }, with the two closure equations: 

V 2 X C V = 2{u + a) C v , 

V 3 X C U = 4u x C u + 2v x C v + 8(u + a) V X C U + 6vV x C v . (44) 



The presence of a in the gauge-invariant coefficients of ( [4"4| ) proves that the 
parameter a is essential, not removable from the studied ZCR by gauge 
transformations. 



10 



Next we pose the problem of finding the hierarchy of all systems 

1H f tjUjVj . . . j U x ...xi V x ...x) i 

v t = g (x, t,u,v,..., u x ... x , v x ... x ) (45) 

which admit ZCRs ([!]) with X determined by ( f43| ) and T decomposable over 
the cyclic basis as 

T = Cl C u + c 2 C v + c 3 V x C u + c 4 V x C v + c 5 V 2 x C u , (46) 

where ci, . . . , C5 are any scalar functions of x, t, u, v, . . . , u x „, x , v x „, x and a. 
Substituting T ( f4"B"D into the characteristic form of ZCRs 

fC u + gC v = V x T (47) 

and using the closure equations (0), we obtain 

ci = -D x c 3 - 8 (u + a) c 5 , 

c 2 = -D x c A - 6fc 5 , c 3 = -D x c 5 , 

f = D x d +iu x c 5 , g = D x c 2 + 2 (u + a) c 4 + 2v x c 5 , (48) 

where the functions C4 and C5 remain arbitrary. The relations (§^) lead to 
the following expression for the right-hand sides of the represented equations 



f ^=(M + \N)p, (49) 

where A = — 8cr, p = (c 5 , — c 4 ) T (T denotes transposing), 

- 8uD x - Au x N 
-QvD x -Av x Dl-2u, 



N=( D Q X iV (51) 



Finally, we impose the conditions df/dX = and dg/dX = on (f49D , in 
order to find the set of admissible two-component functions p. Following the 
same way of reasoning as used in Section [|, we conclude that the systems 
fl4"5] ) represented by ([!]) with X (|3|) constitute a discrete hierarchy generated 
by the recursion operator R = MN^ 1 with M (|50| ) and N ([51]). Inverting 
the matrix differential operator N, 

N~ l =( D ( °X (52) 
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we obtain the explicit expression for the recursion operator R, 

'.Dl - 8u - Au x D- 1 \ 
R=\ \ , , 53 



which coincides with the one found in [11 . 



6 Coupled KdV— mKdV equations 



Recently, Kersten and Krasil'shchik |L3[ introduced the new system of cou- 
pled KdV-mKdV equations 

Ut = —u xxx + 6uu x — 3ww xxx — 3w x w xx + 3u x w 2 + 6uww x , 

w t = —w xxx + 3w 2 w x + 3uw x + 3u x w (54) 

and found its recursion operator (in the sense of higher symmetries). The 



recursion operator of (|54T) has an unusual structure: its nonlocal part consists 
of terms of the form t\D~ x , e 2 D~ 1 e 3 and e 4i D~ 1 e 5 D~ 1 , and those expressions 
ti are polynomials of u, u x , w, w x , w xx , w xxx , y, sin (2z) and cos (2z), where 
y and z are nonlocal variables, 

y : y x = u, z : z x = w. (55) 

The appearance of trigonometric functions of a nonlocal variable in the 
recursion operator derived in [|13[ looks strange and, of course, may be ex- 
plained in various ways. One possible explanation can be obtained by the 
method described in Sections | and [|, which allows us to re-derive the recur- 
sion operator of (|5lf ) in the quotient form R = MN -1 , where M and N are 
2x2 matrix linear differential operators of orders four and two, respectively, 
whose coefficients contain only local variables. This approach makes use of 
the matrix X of the ZCR of (Q), found very recently in ||14j| , 




u - 

X = I 1 a J . (56) 

\0 6m 

where a is a parameter. 

Starting from the characteristic matrices C u = dXj du and C v = dXj dv 
and repeatedly applying the operator V ' x = D x — [X, ], we find the cyclic 
basis to be eight- dimensional, 

{C u , C v , V X C U , W X C V , W X C U , V^Ct,, W X C U , V^C„} , (57) 



12 



for all nonzero values of a. Omitting here the cumbersome closure equations, 
i.e. decompositions of V X C U and V X C V over (|57|) , we only note that their 
gauge- invariant coefficients contain the parameter a, which is an essential 
parameter therefore. Then, the problem of finding the hierarchy of coupled 
evolution equations {u t = f, w t = h} represented by (|I|) with X (|56)) leads 
us to 

^ - (M + XN) ( C7 ) , (58) 



hj \c 8 

where A = 36a 2 ; the unknown functions C7 and c 8 of A, x, t, u, w, . . . , u x _ x 
V0x...x are the coefficients at V X C U and V^C„, respectively, in the decomposi- 



tion of T over (57); and the cumbersome differential operators M and N are 
given below, in a simplified form. Finally, the conditions df/dX = dh/dX = 
lead us to the conclusion that the represented hierarchy is generated by the 
recursion operator R = MN^ 1 . 

The second-order differential operator N can be factorized as 

0\ „ (P xx P 12 



N - DP - °H oj' £>• (59) 



where 



Pa = Dm + «. + ^±2^ > (60) 



'X 

W W x 



u x 

D i2 = -2wD x - 12w x 

w 

.2 1 «„,„., 2 1 c„,,4 



u + Quw z + 8u> 4 - 2ww xx uw xx . , 

+ 77— > ( 61 ) 



w 2 



P21 = — , (62) 

w x 



„ „ uw + 4w; 3 — w xx , . 

P22 = D X ^. 63 

w x 

The fourth-order differential operator M can be expressed as 

M = R loc N + Q, (64) 
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where R\ oc turns out to be exactly the local part of the recursion operator R 
found in (so we are on the right way), 

, -Dl + Au + w 2 -2wDl - w x D x + 3uw - 2w xx . 
^ioc = 9 , , 65 

1 2w -D 2 x + 2u + w 2 1 

and the components of the first-order differential operator Q, 

Q=( Q n 11 % 2 ), (66) 



K Q21 Q22 

are determined by 

Qn = (2u x + ww x ) D x + uw 2 + w 2 x + ww xx 

2u x w x (2uw + 5w ) 
H H , (67) 

w w v 



Q12 = (~3wu x + (2u - Aw 2 ) w x - w xxx ) D x 



u w 



2u 2 

J ■■■ Auw 3 - 23u x w x - 8wwl + (2u - Aw 2 ) w xx 

w 

{2u 2 + 13uw 2 + 20w 4 — 5ww xx ) u x — {uw + Aw 3 — w xx ) w xxx 

w x 

2wu x w xx . ^ 

WW r ' 



2 2 
W W U 

Q21 = w x D x + 2uw + 5w 3 + + (69) 



Q22 = (u x + ww x ) D x — 2u 2 — Y3uw 2 — 20w 4 — 14-u^ + 5ww x 
u x w x — uw xx (uw + Aw 3 — w xx ) u x 



w w 7 



(70) 



Now, let us proceed to the nonlocal part QN 1 of the recursion operator 
R = MN~ l . Since N' 1 = P^D' 1 due to (||), we have to invert the 
operator P determined by (|6"0])-(B3|). Considering 



s = ^' = it Z) : sp = ps = (I ! )• < 71 > 
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we obtain 



Sn — —P21P22S21, 



521 — — {P11P21 1 P22 — P12) , 

522 = {P11P21 P22 — P12) PnPii 1 - (72) 



Then the relations (|72|), (|60|)-(|63|) and (|67|)-(|70D lead us to a cumbersome 
explicit expression for the nonlocal part QP^D^ 1 of the recursion operator 
R in terms of local variables and operators D x , D~ l and K~ x , where the 
operator K~ l is formally inverse to the second-order differential operator 



K — Pi\P 2 i P22 — P12 



w x ^9 /2w„ wz \ „ / w xxx w x w xx \ ,„„. 

= -^D 2 x + _f£-_f }D X + U Wx + -?Zl--Z-£l . 73 

Note, however, that no nonlocal variables have appeared as yet. The 
nonlocal variable z (|55| ) appears only when we make the factorization of the 
operator K fl7£~ 



-2iz „4iz 

e — -2iz 



if = ^ D x e~ 2lz w x , (74) 

in order to express the inverse operator K~ l in a more conventional form, 
i.e. in terms of D' 1 : 

K- 1 = —D-^-^wD-^w. (75) 



w 



x 



Eventually, when we rewrite the nonlocal part of R in the form given in 
IB] , these exponentials e 2lz and e _4lz from K^ 1 fl75] ) give rise to sin (2z) and 



cos (2z), while the nonlocal variable y (|55| ) appears simply as a result of the 
identity u + yD x — D x y = used. 



7 Conclusion 

Throughout the paper, we only used the cyclic bases of ZCRs as a conve- 
nient tool, well adapted to the gauge- covariant nature of ZCRs. We believe, 
however, that the cyclic bases of ZCRs can admit some interesting and useful 
geometric interpretation. 

All the examples, studied in this paper, concerned evolution equations 
and systems of evolution equations. Of course, this is not a sign of a serious 
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limitation for the described technique. It is well known that a very wide class 
of systems of PDEs, including at least all normal systems []15| , can be brought 
into the evolutionary form (see |T^] for a nice explanation of this). There 



are many advantages of using the evolutionary form of a studied system. 
One of them, related to ZCRs, consists in that the cyclic basis of a ZCR 
is constructed of the characteristic matrices and their successive covariant 
x-derivatives only, because the covariant t-derivatives of the characteristic 
matrices (i.e. V t C u = D t C u — [T,C U ], etc.) can be expressed through the 
covariant x-derivatives in the evolutionary case ||]. 

Further examples of using cyclic bases of ZCRs can be seen in [171 ITS, IIP 



In ||17|| , where a ZCR with a parameter was found for the Bakirov system 
that admits only one non-Lie local generalized symmetry, a gauge-invariant 
coefficient of a closure equation was used as an indicator of non-removability 
of the parameter, like we did it in Sections [|-|5] of the present paper. Note, 
however, that other types of gauge invariants can also be helpful in problems 
of gauge non-equivalence of ZCRs and non-removability of parameters ||| . 



In [18], a recursion operator for a new integrable system of coupled KdV 



equations was derived from the eight-dimensional cyclic basis of the sys- 



tem's ZCR. The recursion operator, obtained in [|T^], has a strange non- 
local part, expressed in terms of local variables and operators D^ 1 and 
(3-D^ — 4t> D x — 2v x ) . Note, however, that the nonlocal variable w : v = 
3w xx /w and the factorization 3D^ — AvD x — 2v x = 3w~ 2 D x w 2 D x w 2 D x w~ 2 
make it possible to express this unusual operator (3D^ — ^vD x — 2v x ) 1 in a 
more conventional form, i.e. in terms of D' 1 , like we did it by (|75|) for the 
Kersten-Krasil'shchik recursion operator. 



In [|]j| , the technique of cyclic bases was applied to the problem of how 
to distinguish the fake Lax pairs (introduced by Calogero and Nucci) from 
the true ones. A clear difference was found in the structure of cyclic bases, 
namely, the closure equations in the 'fake' case turned out to be first-order 
and typical for conservation laws, like (|2B|). Though the fake Lax pairs are 



simply gauge-transformed matrix conservation laws, it still remains not so 
clear what are those true Lax pairs which are generally expected to represent 
only integrable (in some reasonable sense) equations. See, in this connection, 
one irreducible sl2-valued ZCR with an essential parameter, quoted in || as 
Example 6. 
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